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ABSTRACT 

The design and development of multispectral remote sensor 
systems and associated information extraction techniques should 
be optimized under the physical and economic constraints encoun- 
tered and yet be effective over a wide range of scene and environ- 
mental conditions. Direct measurement of the full range of con- 
ditions to be encountered can be difficult, time consuming, and 
costly. Simulation of multispectral data by modeling scene, atmos- 
phere, sensor, and data classifier characteristics is set forth as 
a viable alternative, particularly when coupled with limited sets 
of empirical measurements. A multispectral system modeling capa- 
bility is described. Use of the model is illustrated for several 
applications -- interpretation of remotely sensed data from agri- 
cultural and forest scenes, evaluating atmospheric effects in 
Landsat data, examining system design and operational configura- 
tion, and development of information extraction techniques. 


INTRODUCTION 

Remotely sensed, multispectral scanner (MSS) data have been shown to be 
of value in a variety of applications related to earth resources surveys. The 
design and development of future scanner systems and related information extrac- 
tion algorithms must be optimized under a variety of both physical and economic 
constraints to yield maximum amounts of information at a minimum cost for a 
particular application or a series of applications. In addition, optimal infor- 
mation extraction algorithms must be developed to utilize present scanner sys- 
tems (e.g., Landsat) for applications of current interest. Such optimal develop- 
ment is extremely difficult, however, due to the wide range of parameters availa- 
ble to the scanner system designer, the wide range of scene characteristics of 
interest in a particular application, and the wide range of atmospheric condi- 
tions which may be encountered in an operational environment. 

APPROACH 

A viable approach to these difficult design and development problems is to 
simulate MSS data based on theoretical models of scene, atmospheric, scanner 
system, and data processor characteristics. Employing realistic measurements 
of the parameters required by these models, data can be generated at a rela- 
tively small expense over a wide range of conditions of interest. These data 
can then be analyzed for the purposes mentioned above. Simulated data can also 
be used to supplement empirical data obtained by an experimental or operational 
sensor system, to facilitate a greater understanding of those data. 
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A Multispectral System Simulation Model has been developed at ERIM, based 
on previously developed capabilities for mathematically modeling the physical 
processes involved in the generation of remotely sensed signals. Outputs from 
the model have been employed for a variety of applications, some of which are 
illustrated later in the paper. 


DESCRIPTION OF MODEL 

The present basic components of the ERIM Multispectral System Simulation 
Model are; 

(a) bidirectional reflectance model, 

(b) atmospheric radiative transfer model, 

(c) sensor system model, and 

(d) data processor model. 

An overall block diagram of the model is presented in Figure 1. The basic 
equation used to compute the spectral radiance L^ observed by the scanner is: 

= L^(S)t(X) + L^(P) 


where X is wavelength. 


t(A) is the spectral transmittance of the atmosphere along the path 
between surface and sensor. 


L, (P) is the spectral path radiance due to multiple scattering by the 
atmosphere along the path, and 

L, (S) is the spectral radiance emanating from the surface, i.e., the 
radiance before any atmospheric effects occur. 


The surface radiance, in turn, depends on the direction and magnitude of 
irradiance at the surface, the scene bidirectional and diffuse reflectances, 
and the direction of observation. Specifically, the surface spectral radiance 
is : 


L^(S, e^, 0g, i()g) = - E^ (Direct, 6^, • p (Bidirect , X, 6^, 9^, Hj^) 

+ E^ (Diffuse) • p (Diffuse , X) 


where E^ (Direct) 


E^ (Diffuse) 


is the direct (solar) spectral irradiance, 

is the diffuse or hemispherical (sky) spectral irradiance. 


p (Bidirect , X) is the bidirectional spectral reflectance of the surface, 
relative to that of a horizontal, perfectly Lambertian 
surface , 

p (Diffuse, X) is the Lambertian spectral reflectance of the surface 
(i.e., response to a diffuse radiation source). 


and (®v’ "^v^ ^®s’ zenith and azimuth angles of view and sun, 

respectively. 
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Both the irradiance terms and the terms Lp and T depend on the properties of 
the atmosphere for the condition under consideration. In addition, 

and Lp depend on the albedo Pg of the background surrounding the surface element 
being viewed. 

The bidirectional reflectance model used is that developed by Dr. G. H. Suits 
of ERIM [1,2]. This model computes the bidirectional reflectance of a vegeta- 
tive canopy as a function of the spectral reflectance and transmittance charac- 
teristics of each of several canopy components and their sizes and average 
orientations within the canopy (i.e., horizontal and vertical projections of 
their cross-sectional areas) . The basic model assumes that the canopy is infi- 
nite in extent and that its components are randomly oriented in layers. Com- 
puted values have agreed well with measured values in several instances [3-7]. 
Exceptions have been noted for row crops by two investigators [5,6]. Verhoef 
and Bunnik incorporated a crop row-effect parameter in a one-layer Suits' model 
and reported improved agreement between subsequently calculated and measured 
values [6]. In both References 6 and 7, certain discrepancies between modeled 
and measured reflectances have been attributed to specular leaf reflectances. 
Suits believes they could be a consequence of his use of only vertical and 
horizontal elements to represent leaf distributions and is developing a simple 
modification of his model to overcome this problem. Recently, Suits has added 
slope and aspect (i.e., azimuthal direction of slope) parameters to enable pre- 
diction of reflectances from canopies on sloping terrain (an example is pre- 
sented later) . 

In the event that empirical measurement data are available, such as labora- 
tory or field-measured reflectance spectra, these spectra can be utilized in 
simulation calculations in place of spectra computed theoretically with the 
reflectance model. Similarly ,' measured radiance spectra can be inserted at the 
appropriate place in the chain of calculations. 

The atmospheric radiative transfer model is that developed by Dr. R. E. 

Turner of ERIM [8,91. This model effectively has two parts, one which computes 
direct and diffuse solar irradiance spectra on a surface (at sea level or at 
any specified altitude) , and the other which computes the spectral path radiance 
and spectral transmittance of the atmosphere along the path between the surface 
of interest and the sensor. These computed spectral quantities all depend on 
the state of the atmosphere and some depend also on sun position, view geometry, 
and/or background surface, as noted earlier. The diffuse irradiance term com- 
puted by the model includes all radiation scattered at least once by the atmos- 
phere. Usually, this radiation has been assumed to be uniformly distributed 
throughout the hemisphere; however, it would be more accurate to take into 
account its actual spatial distribution which is concentrated in the aureole 
which surrounds the sun. The model assumes multiple scattering in a plane- 
parallel atmosphere. Absorption by aerosols can lie included in the computation, 
although it usually has been omitted in calculations made. The most direct way 
to describe an atmospheric state for the model calculations is to input measured 
atmospheric optical thickness readings. Alternatively, one can utilize standard- 
ized atmospheric states which are labeled by their horizontal visual ranges. 

The major component of the sensor system model is the set of functions that 
represent the relative spectral responses of the scanner's spectral channels. 

In the model, these can have arbitrary shapes as a function of wavelength. A 
radiance responsivity calibration option exists to convert effective inband 
radiances to sensor output units, and a sensor noise option is planned. For 
studies of satellite system operational considerations, an orbital model and 
sun position model permit one to predict system responses for a variety of 
orbit inclinations, equatorial crossing times, and scanning configurations. 

Multispectral data classifier perfoinnance is simulated through a Monte Carlo 
technique based on statistics determined for the spectral classes of interest. 
Random samples are drawn from multiple scene-class populations and classified 
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into specified recognition classes using stati -i'cal decision rules. In addi- 
tion to pure single-class pixels, multiclass mixture pixels which include field 
boundaries and/or spatial misregistration can be simulated [101. 


APPLICATIONS OF MODEL 

The multispectral system simulation capability described above has proven 
to be useful for a number of applications. Several representative uses are 
illustrated and described in this section, beginning with agricultural appli- 
cations . 

Knowledge of the effects of various factors on crop reflectance and their 
relative importance are both important in the design and development of pro- 
cessing techniques and procedures for crop inventory. Winter wheat is a crop 
that has a rather distinctive temporal development pattern relative to other 
crops and natural vegetation in the U.S. Great Plains. Wheat reflectance as a 
function of time was simulated using the model and field-measured wheat canopy 
component characteristics [11]. Seven growth stages, from emergence through 
harvest, were simulated (see Table I). 


TABLE I. WHEAT GROWTH STAGES SIMULATED 


Stage 

Number 

Simulated 

Date 

Wheat Growth 
Stage 

Percent Growth Cov 
(Medium Density) 

1 

Mid November 

Emergent 

14 

2 

Mid April 

Jointing 

44 

3 

Mid May 

Pre-heading 

79 

4 

End May 

Post-heading 

82 

5 

Early June 

Senescing 

64 

6 

Late June 

Ripe 

60 

7 

Early July 

Harvested 

14 


The sun position was chosen to match that at Landsat overpass time on each date 
at 38°N latitude. Atmospheric effects and Landsat inband radiances then were 
computed for several different atmospheric states. The resulting time tracks of 
wheat radiance in Landsat Bands 5 and 6 are presented in Figure 2 for one atmos- 
pheric state. 

The family of plots in Figure 2 illustrates that wheat reflectance depends 
substantially on two canopy parameters -- density of vegetation, which differs 
from column to column in the array of plots, and soil brightness, which differs 
from row to row. The sparse, medium, and dense labels of canopy density apply 
on a relative basis at any given stage of development -- on an absolute basis, 
the dense canopy at emergence or post harvest had a lower percentage ground 
cover than did the sparse canopy at the time of full canopy development (i.e., 
at heading) . Shadowing of exposed soil by the wheat leaves is another factor 
affecting the canopy reflectance. As one would expect, the soil brightness 
affected the total radiance most when ground cover percentage was lowest -- 
these and other simulated values serve to quantify the magnitude of the effect 
and indicate the variety of signature time tracks that are possible. 


A related parametric study examined in more atail the characteristics of 
wheat reflectance during the "greening up" stags, of wheat development [12]. 
Again, Landsat inband values were computed and analyzed. More levels of canopy 
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density were included through a parameter GLAI (Green Leaf Area Index) which is 
the total single-sided surface area of green leaves per unit area of ground. 
Figure 3 presents a family of curves for which GLAI varies from 0 (bare soil) 
to 26 (an extremely dense value to serve as a limiting case) with the mean soil 
brightness; typically mature wheat GLAIs are on the order of two to four, occa- 
sionally as high as six or eight. If GLAI were fixed, the percent ground cover 
and canopy reflectance could both vary, depending on the structure of the canopy 
Canopy structure in the Suits model is implicit in the ratio of vertical to 
horizontal leaf area projections (V/H ratio) which also can be thought of as 
defining an average effective leaf orientation. The effect of varied V/H ratio 
is also displayed in Figure 3. (The actual V/H ratio of a wheat field will 
change throughout its development cycle.) 

Families of curves for light, medium, and dark soil brightnesses are super- 
imposed in Figure 4. For sparse canopies, the soil brightness is the predomi- 
nant factor, while for dense canopies the V/H ratio predominates. Relatively 
small changes in reflectance occur when GLAI increases beyond a value of two. 
Alternative displays of the data in Figure 4 suggested that a polar coordinate 
description of the data would tend to decouple the effects of soil brightness 
and canopy density. Specifically, in Landsat reflectance space, an increase in 
soil brightness causes an increase in the radial component, while an increase 
in GLAI can be correlated with a counterclockwise rotation of a radial line 
through the origin. This rotation is both strikingly independent of soil bright 
ness and not too sensitive to V/H changes. For Landsat signal space, a similar 
result is predicted but with a displacement of the origin of rotation due to 
path radiance effects. The polar angle for Figure 4 would be the angle whose 
tangent is the ratio of Band 7 to Band 5 reflectance, a quantity that has 
been found useful in Landsat analyses as a green vegetation indicator [e.g., 
13-15] . 

The results presented here explain the reasons for the success of the 7/5 
ratio and derivatives thereof and suggest that the polar angle, or a "green 
angle" measured from the radial line of bare soil, might provide a more linear 
measure of green vegetation density. Rather than utilizing just two of the 
four Landsat bands for defining the angular and radial components, linear com- 
binations of all four such as those of the tasseled cap transformation [16] are 
recommended. The green angle is most closely related to the green component of 
the tasseled cap transformation. 

The development of the tasseled cap transformation was based on calculations 
made using the Suits model, as well as analysis of empirical Landsat data. The 
majority of variance in Landsat data lies in the plane defined by the brightness 
and green axes. This transformation has proven to be a useful data processing 
technique and interpretation aid. 

The above calculated reflectance values compared favorably with empirical 
values measured with a Landsat-band radiometer in wheat fields throughout a 
Kansas growing season, as shown in Figure 5. 

Forests are frequently found in mountainous or hilly regions on sloping 
terrain. Slope and aspect do have pronounced effects on remotely sensed data. 
The slope and aspect parameters of the Suits model were used to predict the 
reflectance of a pine forest under various slope and aspect conditions. 

Figure 6 presents one set of results for Landsat Band 6 which indicated the 
extent to which the reflectance of a pine canopy is predicted to depart from 
that of a Lambertian (perfectly diffuse) surface. 

Atmospheric effects are present in remotely sensed data acquired for all 
applications. An understanding of the magnitude and character of these effects 
is important to both the interpretation of data and the development of informa- 
tion extraction techniques and procedures. Figure 7 (reproduced from Reference 
17) illustrates several facts about Landsat Band 4 data. First, path radiance 
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constitutes a large fraction of the signal from a target of medium reflectance 
('''8%). Second, the amount of atmospheric haze (indicated by horizontal visual 
range) affects the signal received. Finally, there are appreciable scan-angle- 
dependent effects even over the limited ±6° scan range of the Landsat MSS. Cal- 
culations of atmospheric effects can be tailored to the particular conditions 
and parameter values appropriate for any given application study. 

The system model also has been used to examine the influence of scan geome- 
try and system orbit parameters and sensor spectral characteristics on scanner 
responses [18]. In addition to atmospheric effects, the question of surface 
glint effects was addressed. Figure 8 presents a map of orbital coverage from 
70°N to 70°S latitude along a ground track simulating the Thematic Mapper on an 
earlier proposed orbit with an 11:00 AM equator crossing. The cross-hatched 
zones indicate the surface slopes for which direct specular glint conditions 
could exist (e.g., from the sea surface) as a function of satellite latitude. 

By way of explanation, surface areas in the zone labeled >.15° would be suscepti- 
ble to glint from any surface slope of 15° or greater; the enclosed zones for 
.>10° and >5° would be included as well in the potential glint region. The pre- 
dicted glint potential for the currently planned 9:30 AM Thematic Mapper orbit 
is much less, with surface slopes of at least 20° being required and the >.20° 
zone roughly matching the >.10° zone of Figure 8. 

Figure 9 presents the nominal spectral bands of Thematic Mapper and Landsat 
MSS along with a reflectance spectrum of a green leaf to indicate the potential 
utility of the model in addressing problems of sensor spectral band selection. 
These problems require a consideration of atmospheric and viewing conditions 
and effects in addition to the spectral characteristics of scene surfaces. 

Data processor simulation, based on signature statistics, is exemplified by 
previously reported predictions and analysis of the effects of spatial misregis- 
tration on classification performance [9]. 


SUIIMARY 

The ERIM Multispectral System Simulation Model has been described and 
several illustrative applications have been presented which address agricul- 
tural, forestry, general Landsat data interpretation and analysis, and system 
parameter selection problems. The capabilities described represent one imple- 
mentation of a class of physical models which can play an important role in the 
continued development of remote sensing technology, especially when coupled with 
appropriate empirical measurements . The full potential of such modeling is yet 
to be realized. 
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FIGURE 1 


FIGURE 2 


ERIM MULTISPECTRAL SYSTEM SIMULATION MODEL 



TEMPORAL TRACKS OF SIMULATED WHEAT CANOPY 
RADIANCE SIGNATURES IN LANDSAT BANDS 
(Key to growth stages is in Table I). 
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FIGURE 7 


FIGURE 6 

Interpretation of Data 
From Forested Areas 
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Prediction of Atmospheric Effects 
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FIGURE 8 

Sensor Orbit Selection 



FIGURE 9 

Sensor Spectral Band Selection 
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